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An accurate measurement of the bunching of photons in the fluorescent emission 
from an ultracold ensemble of thermal ^'^Rb atoms in a steady-state magneto-optical 
trap is presented. Time-delayed-intensity-interferometry (TDII) performed with a 5- 
nanosecond time resolution yielded a second-order intensity correlation function that 
has the ideal value of 2 at zero delay, and that shows coherent Rabi oscillations of 
upto 5 full periods - much longer than the spontaneous emission lifetime of the excited 
state of Rb. The oscillations are damped out by ~ 150ns, and thereafter, as expected 
from a thermal source, an exponential decay is observed, enabling the determination 
of the temperature of the atomic ensemble. Values so obtained compare well with 


those determined by standard techniques. TDII thus enables a quantitative study of 
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the coherent and incoherent dynamics, even of a large thermal ensemble of atomic 


emitters. 


PACS numbers: 
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Intensity interferometry, or the measurement of photon correlation, provides a wealth of 
information regarding the light source and the mechanism of emission. A study of the mere 
arrival times of photons from the source enables one to determine its nature - whether it is 
thermal (chaotic), or coherent, or quantum. The form of approach to the asymptotic value 
allows one to infer further details of the source. For example, one of the earliest measure¬ 
ments [T] of spatial photon correlation of light from the star Sirius enabled the determination 
of the diameter of the star. Intensity interferometry is routinely used in particle physics PE] 
to study decay processes and to deduce the interaction between particles. Other applications 
include the search for naturally occuring non-classical sources of radiation in astrophysics 
[1] , study of light emission from nanostructures [5] and particle size measurements [6] . Typ¬ 
ically, pairs of detectors either seek the simultaneous arrival of photons, or measure the 
delayed arrival of a photon at one detector with respect to the arrival of a photon at the 
other. The measurements are quantihed by the second order correlation function, (also 


known as intensity-intensity correlation function), r) = 


+ R,t + r)^ 

(^/i(r,f))(/2(r,t)^ 

where and Rir + RR + t) are the intensities of light reaching detectors D1 and D2 at 


locations r and r-|-R at times t and t + r and the angular brackets denote time averaging. It 


is the correlation of the number of photons, or intensities, that is examined, as opposed to a 
correlation of amplitudes in conventional interferometers, and hence the name time-delayed 


intensity-interferometry (TDII). It has been theoretically shown that (see, for example |2]) 
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that the second order correlation fnnction = 0) = 2 for a thermal state, implying 

a tendency for bunched or correlated emission of photons; is unity for a coherent source, 
implying emission of photons at random times, and equals 1 — ^ for a n-photon Fock state, 
signifying anti-bunching. The value of for all sources, however, approaches unity for 

long time delays. 

Temporal bunching of photons from thermal sources, ever since the postulation of the 
concept, has been an intriguing phenomenon, and has been the focus of numerous experi¬ 
ments. The earliest laboratory thermal source studied was a Hg vapour lamp [^, light from 
which showed a meagre bunching of 1.17. Martiensen and Spiller devised a method of creat¬ 
ing pseudo-thermal light by transmitting coherent laser light through a rotating ground glass 
plate [S] such that the time-varying surface inhomogenieties introduced temporal and spatial 
decoherence. In recent years, laboratory control and measurement techniques have enabled 
creation of pseudo-thermal light sources with theoretically expected values of = 2 

[SHig. 

We report here Time-Delayed Intensity-Interferometric measurements on light from an¬ 
other source of bunched photons - an ensemble of laser cooled atoms. Though laser cooled 
atoms have been available for more than three decades, direct TDII measurements of their 
emission have been very few fTBHT^] . All measurements hitherto have been carried out in op¬ 


tical molasses, which were either periodically, or continuously loaded with atoms precooled 
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in a MOT. In this paper we present measurements of the second order correlation func¬ 
tion of light emitted by ultracold atoms in a steady-state magneto-optical trap, where the 
cooling and repumper beams, and also the quadrupolar magnetic held are kept on. We ob¬ 
serve the ideal value of 2 for the zero-delay intensity-intensity correlation function. Damped 
Rabi oscillations are observed for time delays upto ~ 150ns, and an exponential decay for 
longer time delays. Despite the fact that the emission being studied is from a collection 
of uncorrelated atoms that are in random thermal motion, and that the observations are 
averaged over an 8-hour period, coherent effects are seen, bringing out the power of higher 
order correlations in revealing hidden periodicities and providing a measure of coherent and 
incoherent dynamics. The exponential decay at long time delays was used to determine the 
temperature of the ensemble. 

Time delayed Intensity Interferometry was performed on the huorescent emission from 
®^Rb atoms cooled and trapped in a magneto-optical trap (MOT) (Fig. [^, which differed 
from usual MOTs, in that the two pairs of beams in the x-y plane were steeply inclined to 
each other, enclosing an angle of 55° rather than the usual 90° so as to accomodate, within 
the chamber, a pair of lenses of short-focal length and high numerical aperture. These lenses 
were positioned facing each other, such that their focal points coincided with the centre of 
the MOT and could thus be used to focus light onto the MOT, or collect light emitted by a 


small volume within the cold cloud. To avoid clipping at the lens mounts leading to unde- 
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FIG. 1. Schematic of the experiment showing the collection of cold atoms inside the MOT chamber. 
A lens within the chamber collects fluorescence from the atoms, which is passed onto the TDII setup 
where the arrival of individual photons are time-stamped and archived. Bold arrows show the cooling 
beams in the x-y plane. The z-beams (perpendicular to the page) are not shown. 
sired scattered light, the diameters of the beams in the x-y plane were restricted to ~ 1.5mm 
while the z-beam had a diameter of 8mm. The cooled and trapped atoms were viewed using 
a CCD camera, and their number estimated by collecting part of the fluorescent light onto 
a femtowatt detector. The typical cloud was roughly ellipsoidal with a mean diameter of 
~ 400/im, and contained about 20000 atoms. 


An intensity-interferometer was formed using a hber splitter, where the two output ends 
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were connected to two high speed APDs that had quantum efficiencies of 65% and a dead 
time of 30ns. The output of the APDs were fed to a homebuilt, FPGA-based time-tagged- 
single-photon recorder, APODAS (Avalanche Photodiode Optical Data Acquisition System 
pT] that utilised high speed ethernet connectivity and stored, in a PC in realtime, the arrival 
times of all detected photons with a temporal resolution of 5 nanoseconds. Post-processing 
by software enabled the determination of for all r from a single recording of the 

data As we worked in the photon-counting mode, the expression for in terms of 

coincidences is |T%] 

g^^^T) = N,T/{N,N2T,) ( 1 ) 

where r at is the delay between arrival at the two detectors, A^i, N 2 and Nc are the number 
of counts at detector 1, detector 2 and the coincident counts respectively. T is the total 
observation time and Tc the time window for coincidence (arrival of two photons is considered 
simultaneous if they are detected within a time gap of Tc). 

For determining the bunching characteristics of emission from an ultracold atomic en¬ 
semble, ®^Rb atoms were laser cooled from close to room temperature, to ~ lOO/iK. The 
MOT was extremely stable, with the lasers locked and the cold cloud in steady-state for 
days. A typical run of the experiment lasted 8 hours. The cold cloud was obtained, and the 
cooling and repumper beams, and the magnetic held were kept switched on for the entire 


duration of the experiment. The cold ensemble was constantly monitored by imaging it on 


a camera, and also by measuring the fluorescence on a femtowatt detector (see Fig. [^. For 
the purpose of determining its second order correlation function, fluorescence light from the 
central region of the cloud was collected by the high-numerical aperture lens placed inside 
the sample chamber, a few millimeters from the trap centre. Care was taken to ensure that 
no part of the laser light entered this lens, either directly, or upon being scattered by the 
parts of the MOT chamber. The light was conveyed by a series of mirrors and lenses to 
the input of the TDII setup. The count registered in the presence of the cold cloud was in 
the range 40,000 - 80,000/s while in the absence of the cold cloud it reduced to ~1200/s, 
confirming that it is predominantly light from the cold atoms that enters the TDII setup. 
Photon arrival time data was recorded for various detunings of the cooling beam. For each 
detuning the number of atoms trapped was estimated from the fluorescent intensity recorded 
on the femtowatt detector. While the total number of atoms collected ranged from ~ 8000 
to ~ 22000, the high numerical aperture lens employed for accepting light for TDII measure¬ 
ment restricted the collection of light to that from approximately one-tenth of the volume of 
the cold ensemble. The temperature of the collection of atoms was determined by the trap 
oscillation method m and was found to range from 200/rK to 50/iK for detunings of the 
cooling laser varying from -12MHz to -22MHz. The effective Rabi frequencies ranged from 
25MHz to 40MHz. 


Two-time intensity-correlation values were derived from the TDII measurements. Fig. 
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FIG. 2. (a) Experimentally obtained vs r for laser light measured before its entry to the 

MOT chamber (red) and for fluorescence light from the cold atoms (brown), and (b) Oscillations 
seen g^"^^ (r) at short time scales, for different detunings of the cooling beams; the curves are drawn 
as a guide to the eye, and the values of the detuning are given in the legend. 

shows the second order correlation functions of the input laser light determined before 
it enters the MOT chamber, and of the light emitted by the ultracold atoms trapped in the 
MOT. It is clear that absorption and re-emission by the thermal ensemble of atoms has led 
to bunching in light, which was initally coherent = 1 for all r). Several interesting 

features are observed in for light from the cold atoms. Periodic oscillations, are seen; 


these are damped out by ~ 150ns, and thereafter the curve decays steadily towards unit 
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value. The curve is less noisy for small r than for larger ones. As the cooling beams 

are further detuned from the 5 S 1 / 2 F = 2 —)■ 5 P 3 / 2 F = 3 transition (the so-called ’’cooling 
transition”), the value of approaches closer to 2 while the oscillations in correlation 

are more prominent and more rapid, and their damping is slower (Fig. |^). 

Let us, for simplicity, assume the Rb atom to be a 2- level system. Irradiation of an atom 
by coherent, near resonant light causes two processes to occur. One is the periodic absorp¬ 
tion and coherent emission of radiation leading to an oscillation between the excited and 
the ground state at the Rabi frequency - a rate determined by the intensity and detuning 
of the incoming radiation. The other is the absorption and spontaneous random emission, 
at a rate that falls exponentially with a characteristic lifetime, which, for the excited state 
of ®^Rb is ~28ns. The atoms in the ensemble, being uncorrelated, do not emit in unison, 
and thus no periodicity will be evident in the direct observation of emission from the col¬ 
lection. However, all atoms undergo Rabi oscillations at the same frequency, and thus have 
a high probability of emission at the same regular interval. The second order correlation, 
which is a measure of the probability of emission at time t + t conditioned on an emission 
having occurred at time t, will therefore exhibit periodic maxima at regular intervals tr, 
the inverse of the Rabi frequency. Thus, the two-time intensity correlation measurement is 
a simple yet powerful technique that can reveal hidden periodicities. The periodic oscilla¬ 


tions arising from coherent dynamics under steady-state driving helds, nevertheless, show 
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decay as decoherence sets in due to spontaneous emission and inter-atomic collisions. This 
reasoning also explains why is less noisy for short time delays than for larger delays. 

The Rabi frequency in our experiment being a few tens of MHz, coherent oscillations should 
occur at the time scales of few tens of nanoseconds. The lifetime of the excited state is 
~ 28ns, implying that a decay in amplitude of oscillation will occur over ~ 100ns ( a few 
lifetimes). Interatomic collisions occur at yet longer time scales. Thus, short delays show 
cleaner curves. 

Let us consider a thermal collection of N independent atoms under the action of a near¬ 
resonant driving held of Rabi frequency fl. The electric held T(r) due to coherent emission, 
at a point of observation r is the resultant of contributions from each atom, j, and may be 
written as (Eq. 13.48 in [18]) 

N N 

~ ( 2 ) 

i=i j=i 

where Sj{e g), is the de-excitation operator for the two-level atom, the dynamics of which 
is given by the master equation for the driven two-level atom (Eq. 13.1 in [IH|), and cfyj the 
phase of the electric vector (due to the coherent emission of the atom) at the point 
of observation (detection) depends on the location of the atom and the orientation of the 
atomic dipole. In a MOT, the resultant driving held due to the six cooling beams varies 
in a complex manner in intensity and polarization from position to position. Likewise, the 


orientation of each atom varies as it moves within the MOT region (see for example [T9]). 
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The second order intensity correlation function is defined as : 

_ (£t(0)£t(r)£(r)£(0)> 

® ^ (5t(0)£(0)>2 

Substituting from Eqj^ to Eqj^ yields terms of the following forms, with appropriate pref¬ 
actors : 

N 

(a) J]{£t(0)fj(r)£,(r)f,(0)) 

i=i 

N N 

('') E E (d(o)fi(T)£,(T)5,(o)) 

i=l j=l,j^i 
N N 

(c) E E (d(o)£i(^)^-(^)^j(o)) 

i=l j=l,j^i 
N N 

(d) ^ ^ {5/(0)£t(r)5,(T)i-,(0)> 

i=l j=l,j^i 
N N N 

WEE E (£’/(0)Tj(r)£^fc(r)£’j(0)) and various permutations 

i=l j=l,j^ik=l,k^i,j 
N N N N 

(/) E E E E (d(0)£l(r)£t(r)£,(0)> 

i=l j=l,j^i k=l,k^i,j l=l,l^i,j,k 

Terms of the form (a) represent single-atom contributions. At r = 0 these show antibunching 
- an atom cannot emit more than one photon at a time. 

Recognising that the emitters are a thermal collection of uncorrelated atoms, the operators 
in the terms in (b) may be re-ordered and factorised to yield : 

N N 

('>) = E E l(d(0)£i(0))(£i(^)^^W>l 

i=l j=l,jj^i 
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where I is the total intensity at the detector due to N atoms. On similar lines, terms in (c) 


lead to the auto-correlation 


(c)4iV(]V-l)|{f;(r)£,(0)>|- 


Terms consituting (d) are related to the anomalous correlation, which, for a thermal cloud, 
vanish on time averaging. Similarly, terms (e) and (f), due to the random phases, also drop 


out upon time averaging. From these arguments, we now hnd that 


(£t(0)ft(T)£(T)f(0))=/V|{fj(0)fj(r)f,(T)£,(0)>| + N(N-i)' 


P , /y/V(f;(T)£,(0)>|2 




Denoting by X = / /iV, is the intensity due a single atom, one obtains 


^ (gj(0)£j(T)g,(r)g,(0)> _ yN r irWfMI" 

n) [ 


g" '[T) = 


+ 1-TT 1 + 


The first term represents the single-atom contributions, which diminishes when the number 
of atoms becomes large. In this limit, the above equation leads to the well known relation 


between the hrst and second order correlations 




As is well known from the Wiener-Khinchtine theorem, the fourier transform of g^^\T) yields 
the power spectral density of the emission. In the time domain it leads to Rabi oscillations. 


/^)(r) = Ao + 



14 


Here + 52 ^ 1 / 2 ^ is the effective Rabi frequency at detuning S. The second order 

correlation function, given by Eq. would then have the form 


g(^\T) = Ao + 

+ + incoherent terms 

Thus, the emission from a collection of N atoms will display coherent Rabi oscillations that 
decay at rates indicative of the relaxation mechanisms [T 8 ] . 

Indeed, our TDII measurements of the fluorescent emission from the cold atoms exhibit 
such oscillations (Figj^b)) - for large detnnings, five full oscillations are seen while for 
small detunings barely one or two are. That the oscillations are reduced in prominence as 
detuning decreases may be understood in terms of the temperature of the atoms. Small 
detunings of the cooling beams result in a hotter collection of atoms, and therefore result 
in increased inter-atomic collisions that decohere the system rapidly. It is thus evident 
that TDII measurements can help determine the relative strengths of various relaxation 
mechanisms as functions of different physical parameters. 

For the collection of cold atoms in our experiments, signatures of the coherent processes 
in the die out within delays of ~ 150ns, and spontaneous emission is expected to have 

caused atoms to make transitions to the ground state within a few lifetimes of the excited 
state. For delay times larger than this, is dominated by effects due to the scattering 
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by moving atoms. While the collection of atoms under study is cooled to ~ lOO/iK, where 
the Doppler width reduces below the natural linewidth, it may seem surprising that the 
effect of the velocity distribution is seen in the scattering. Once again, the power of the 
second order correlation becomes evident, as may be interpreted as the measure of 

the probability of detecting a second photon scattered by an atom with velocity v, within 
a time r of having detected one such photon. When the velocity spread of atoms is large, 
as at higher temperatures, the probability for such an event is low, and thus will fall 

more rapidly towards unity compared to the case for lower temperatures. Thus, at 

large time delays (in this case delays larger than ~500ns) may be used to determine the 
temperature of the ensemble. The elastically scattered light has a Doppler profile determined 
by the velocity distribution of atoms [a EQ]. In the six-beam conhguration, denoting by 
aj{= 2(1 — cosOj) ) the dependence of the Doppler spread of the beam on its scattering 
angle 9j, and by Aj the weight factor for the beam appropriate for its intensity and its 
polarisation and angle dependent elastic scattering cross-section [HEQ], 


6 

9 d\'^) = ^Aiexp{-aj 
i=i 


oolksTr^ 

2mc^ 


(7) 


Here 5 '^^('r) represents the Doppler contribution to the first order correlation function and uio, 
c, ks, T and m are the frequency and speed of light, the Boltzman constant, the temperature 


of the ensemble and the mass of the atom, respectively. Using this in conjunction with the 
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relation E] 




1 + S | 5 '^^(t)P for large r 


( 8 ) 

(9) 


where S depends on the spatial coherence of the light detection system, the temperature 
T may be estimated from the experimentally obtained time-delayed intensity correlation 
function. Fig. [^displays the experimentally obtained valued for as function of r for 

different detunings 6 of the cooling beam. The temperature of the ensemble is determined for 
each detuning of the cooling beam by fitting the experimental data with the corresponding 
(r) curve obtained using Eq. and Eqg As seen in the figure, the data for the different 
detunings £t quite well with the respective curves. It may be noted that the same parameters 
(S', Q!j,Aj) are used for all curves. The values of temperature thus obtained (Ttdii), on 
comparison with the temperature obtained by the trap oscillation method {Tto) show fairly 
good agreement (Figj^, with Ttdii being slightly lower in all cases. 

We now turn our attention to the remaining observation - the increase in with 

detuning. We attribute this to the timing resolution in our experiment. The larger the 
detuning, the colder the collection of atoms, and hence the slower the decay in coherence. 
Thus the timing resolution of 5ns appears adequate for large detunings. For hotter atoms 
obtained at lower detunings, the averaging effect of the time bin becomes discernible, as it 


^ For finite (non-zero) size of source and detector, a factor S is introduced in Eq. ^ (see, for example, m ) 



6 

'1 


4 
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FIG. 3. vs T for different detunings, 6 (1 = -14-5MHz, 2 = -16.1MHz, 3 = -18MHz, 4 = 

-19.4MHz, 5 = -22MHz) of the eooling beam from the cooling transition. The values obtained from 
TDH measurements are shown as dots. The solid curves represent the values of g^"^^ (r) obtained 
from an d\E with temperatures WOfvK, 80fj.K, QOfrK, 50fiK, A3p.K, for the 5 values of detun¬ 

ings, (1 = -14.5MHz, 2 = -16.1MHz, 3 = -18MHz, 4 = -19.4MHz, 5 = -22MHz) respectively. The 
data for all sets are fit with the same choice of parameters, to obtain the temperature. 
is now a larger fraction of the coherence time. 

As pointed out earlier, though light form a thermal collection of atoms is expected to 
exhibit a second-order intensity correlation of 2, this value could not be experimentally 
obtained till very recently; in fact, the present work and that of Nakayama and coworkers 
[15] are the only two reports of this. Further, there has been skepticism on being able to 
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FIG. 4. Comparison of the temperature of the cold atoms, obtained from two methods - TDII 
(circles) and trap oscillation (squares) - as function of the (red) detuning of the cooling beam. 

obtain a good measure of bunching, and of being able to see coherent effects in TDII mea¬ 
surements from a collection of large number of atoms [19]. Several factors have contributed 
to the difficulty in observing the theoretically predicted behaviour. All researchers have 
stressed on the need for good timing resolution, which in the present case is 5ns. However, 
small time bins necessitate longer acquistion times to obtain good statistics, making the 
experiment long and tedious. Light has to be collected over a single spatial coherence region, 
contributing to further reduction in photon counts. These factors require the experimental 
setup to be extremely stable, and the conditions repeatable over nearly ten hours. Another 
factor known to degrade the observation of bunching is the presence of a magnetic field d], 


because of which all measurements hitherto had atoms cooled in a MOT and then trans- 
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ferred to a molasses, either by switching off the magnetic held ([lani]), or by transporting 
the atoms to another vacuum chamber [T3]. Our experiment, however, has been carried out 
in-situ, with atoms in a MOT, with all cooling and repumper beams and the quadrupolar 
magnetic held present. Good mechanical isolation of the setup and temperature stability 
of the environment ensured that the lasers remained locked for the entire duration of the 
experiment. Constant monitoring of the MOT huorescence allowed for corrective measures, 
which, however, were not required. The dihraction-limited collection lens placed within the 
MOT chamber, in close proximity to the cold atoms, and the subsequent spatial hltering 
enabled us to collect light from a small region of the MOT, over which the magnetic held 
was uniform within 2mG, eliminating broadening due to Zeeman shift. Likewise, the low 
temperature, and the thus the sub-natural Doppler width ensured that the Rabi frequency 
is the same for all atoms. Further, the small size of the cold cloud (400/xm across), the 
low number of atoms ensured that reabsorption of the emitted light was negligible. This 
allowed us to detect coherent ehects like Rabi oscillations. In an earlier study, single atom 
dynamics was probed by photon-photon correlation in an optical dipole trap [19], where one, 
two, or three atoms were held trapped. In the present experiment, the number of atoms 
contributing to the collected light is three orders of magnitude higher. Further, atoms move 
in and out of the region from which light is collected, due to the thermal motion, and the 


superimposed trap oscillation in the quadrupolar magnetic held. The transit time of an 
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atom (in the absence of a collision that expels it from this region), is estimated to be ~ lOps. 
The power of TDII is brought out in the present study, where, despite the sample being a 
thermal collection of several thousand atoms, coherent dynamics are revealed. 

In conclusion, we have perfomed Time-Delayed Intensity Interferometry with light emit¬ 
ted by an ultracold atomic ensemble in a steady state MOT. The collection of cold atoms 
is a source of bunched light, where bunching is introduced by spontaneous emission. Well 
dehned, but decaying Rabi oscillatons were seen at small time delays (< 150ns) that give 
way to an exponential decay at larger time delays. It is thus seen that TDII measurements 
enable the study of coherent and incoherent dynamics of the system, providing a relative 
measure of the various dynamical processes occuring at different time scales, even from a 
thermal ensemble of a large number of independent atoms. 
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